Quantization of Spin Direction for Solitary Waves in a Uniform Magnetic Field 
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It is known that there are nonlinear wave equations with localized solitary wave solutions. Some of 
these solitary waves are stable (with respect to a small perturbation of initial data) and have nonzero 
spin (nonzero intrinsic angular momentum in the center of momentum frame). In this paper we 
consider vector-valued solitary wave solutions to a nonlinear Klein-Gordon equation and investigate 
the behavior of these spinning solitary waves under the influence of an externally imposed uniform 
magnetic field. We find that the only stationary spinning solitary wave solutions have spin parallel 
or anti-parallel to the magnetic field direction. 
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I. INTRODUCTION 

Stable localized solitary wave solutions of nonlinear 
wave equations are known to exist ([l|, 0], @, [U, @, @, 
Q). In particular, the existence of stable solitary wave 
solutions has been proven for certain nonlinear Klein- 
Gordon and Schrodine cr equations (0, 0, 0, 0) and 
it has also been established that these solutions can have 
nonzero spin (intrinsic angular momentum) (0, 0). 

Here we consider the nonlinear Klein-Gordon equation 
(NLKG) 



u tt - Au = g(u) 



(1) 



N with N even (see [f|) and the 
R N is defined by g(y) = g(\y\)fa, 
— > R being a continuous function 



where u : R 3+1 — > I 
nonlincarity g : M. N — 

for y ^ 0, g : [0, oo) 
and 5(0) = 0. 

In this paper we will examine the Noether conserved 
quantity 



S[u) = Uf (X x Vu)d 3 X. 



(2) 



called spin, which results from the rotational invariancc 
of NLKG. (We note that for a function u : K 3+1 — > R N , 
we define a counterclockwise rotation of u about an 
axis through the origin in R 3 through an angle 9 to be 
u(Rj X,t) where Rg 1 is a 3x3 counterclockwise rotation 
matrix.). This functional, S, gives the angular momen- 
tum about the origin of a solution u. Our goal will be to 
find the spin of stationary solitary waves when exposed 
to an external uniform magnetic field, B. The notation 
S x , S y ,SLiid S z will be used to denote the magnitude of 
the x, y, and z components of S respectively. 

Equation |T]) can be written compactly using relativis- 
tic index notation as 



d a d a u = g(u). 



(3) 



4), and d a 



d_ d_ d_ d_ 

dx ' dy ' dz ' dt ' 

We model the imposition of 



where 

(__d d_ _d_ d_\ 

\ dx' dy' dz' dt) 

an external uniform magnetic field of strength B = \B\ 
parallel to the z-axis by making the minimal-coupling 
substitutions d a i — > d a — aA a and d a i — > d a — oA a 
into (J3]) giving us 



(d a - aA a ){d a - aA a )u = g{u) 



(4) 



where a is a fixed N x N real skew-symmctric matrix 
with a 2 = — / and 



.4 
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(5) 
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It will be assumed throughout that B is small. This will 
allow us to simplify matters by ignoring terms in Q that 
involve B 2 . We are going to look for stationary solitary 
waves solutions of (0J that carry nonzero spin. A class of 
such solitary waves will be examined and it will be shown 
that they have spin either parallel or antiparallel to the 
uniform magnetic field. 

We look for standing-wave solutions of the form 
u(X,t) = e tn v(X) where Q is some N x N skew- 
symmetric matrix that commmutes with a. We are 
interested in solutions that have v(X) = if?(X)w(r), 
where w[0, oo) — ► M, with r = (x 2 + y 2 + z 2 )?, and 
^ : R 3+1 — > M. is a unit-vector-valued eigenfunction 
of the spherical Laplacian. The Laplacian A can be 
decomposed into radial and angular components A = 
Afi + T-Ajsf. The spherical component Ag, is a second- 
order derivative operator with only angular derivatives 
that acts on real-valued functions defined on the unit 
sphere. It also acts on real-valued functions defined 
on R 3 , leaving the radial dependence unchanged. It is 
known that there exist unit-vector-valued cigenfunctions 
of As with any of the eigenvalues w = —1(1 + 1), where 
I = 0,1,2,3,... (@, 1). Given a nonnegative integer I, 
such an eigenfunction $ : S 2 — > S ,w ~ 1 with eigenvalue 
M exists provided that N > 21 + 1. So A s * = -1(1 + 1)-$, 
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where we extend the action of A5 to vector IR^-valucd 
functions by allowing the operator to act componentwise. 
Thus the coordinate functions of ^ are eigenfunctions of 
A5. Substituting the ansatz 



tively. Thus we have: 



u(X,t) = e tu ^(X)w(r) 
into (1) produces the ordinary differential equation 



(6) 



w"{r) + -w'(r) 
r 



1(1 + 1) 



w(r) + f(w(r)) = (7) 



where I = 0, 1, 2, 3, ... and f(y) = g(y) +0J 2 y with appro- 
priate conditions on / to ensure the existence of localized 
radial soutions (see [j|). It is shown in Q that there 
then exist functions u>[0,oo) — ► R that are exponen- 
tially decreasing far from the origin. Consequently, for 
such functions w, u(X,t) = e ^f(X)w(r) is a solution 
to NLKG. 

We will prove the following theorems: 

Theorem 1. Let B be a fixed real number, a an NxN 
real skew- symmetric matrix, and g : R N — ► M. N where 
g(v) = g(\v\ 2 )-^ with g(0) — for some continuous func- 
tion g : [0,oo) -> R. Let u : K 3+1 -> R N be a solution 
of 



u tt - Aw - Ba(xd y - yd x ) = g(u) 



(8) 



of the form u(X, t) = e tn ^>(X)w(r) where a is an N x N 
real skew- symmetric matrix with the property [Q, a] = 0. 
^ is a unit-vector-valued eigenfunction of the spherical 
Laplacian and w is a solution to ([?]) as described above. 
Define the spin functional S by 



S[u] = u t -{X x Vu)d 3 X. 
IfB^O then S x [u] = S y [u] = 0. 



(9) 



Theorem 2. There are solutions with spin parallel 
(anti-parallel) to B. 

In section II we prove Theorems 1 and 2. Simply put, 
Theorem 1 states that solutions with the ansatz (|6]) to 
equation ^ cannot have nonzero spin components that 
are not parallel to the magnetic field unless B = 0, while 
Theorem 2 imples that the only solutions of the ansatz 
to the equation have nonzero spin parallel to the B field. 

It is important to note that the discussion here is not 
a quantum mechanical one. Although many of the con- 
structions have analogues in quantum mechanics, the in- 
terpretations are different. 
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L y = -i 


f d 

y dx 




L z = -i 


(4 


Inserting 


i) into © 


gives 



d_ 

9y 
d_ 

dz 
d_ 

1 dx 



(10) 

(11) 
(12) 



n 2 <f(X) - Ba(xd y - yd x )V(X)\ w(r) - 
A R w(r) - l ( l+ 2 ^ w(r) + g(w(r))) (JC) = (13) 



Proof of Theorem 1: Since both 51 and a are skew- 
symmetric, then they are both normal. By hypothesis, 
[fi, a] = 0. Thus there exists an orthonormal basis of 
C N consisting of vectors which are eigenvectors for both 
fl and a (0). Also since As and L z have a common 
set of orthonormal eigenfunctions, ...,£_i,£o) £ij ---j £i 
([Tot) tnat s P an the space of all eigenfunctions of As with 
eigenvalue —1(1 + 1), there exists aj m £ C, —l<m<l, 
1 < j < N, giving the expansion 



9(X) 



N I 

E E a ? 

j — l m— — l 



i4>j£m(X) 



(14) 



where {<fij} are an orthonormal basis of C N of eigenvec- 
tors of a and f2 with a<pj = tji4>j where 6j = ±1 and 
Q<f)j = ivj4>j with real Vj. Substituting (|14[) into Q13p we 
get 



N I 



E E «im^fm{(-^ + Bejm)tu(r) -r;(r)}= 0(15) 



j — l m— — l 



where r\ = Anw(r) — ^^"^ w(r) + g(w(r)). Therefore, for 
each j and m either 



a jm = 



or 



{—Vj + ejBm)w(r) — ?/(r) = 



(16) 



(17) 



Next we compute the quantity S x [u] + iS y [u\. It 
should be noted that since ut is real, it is equal to 
its complex conjugate. This will enable us to simplify 
things later. From [T(j we recall that (L x + iL y )£ m = 
\Jl(l + 1) — m(m + l)£ m +i. Using this 



II. PROOFS OF THEOREMS 1 AND 2 

In what follows below, we will denote the x, y, and z 
components of (— i)(X x V) by L x , L y and L z respec- 



N l-l 



S x [u] +iS y [u]=Y^ E ( / w2 i r > 2 dr) x 

j=l m=-l ^° 

VKl+l)-p(p+l) (18) 
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If S x [u] + iS y [u] is to be nonzero, then for some m, j 
and p, Qfj'(p_|_i) ^ and otj P ^ 0. From (flTi)) and (jTTJ) then 



(-i^ + e i B(p+l))u;(r)-T?(r) = (19) 



and 



[—Vj + £jBpjw(r) — i](r) = 
Subtract to get 

Bw(r) = 
=>• B = or to(r) = 

So if S ^ then S x + iS y = 



(20) 



(21) 
(22) 



Proof of Theorem 2: Without loss of generality, we 
assume B > and take 



0-100 
10 
-1 
10 

It is straightforward to show that 



*P0 = - . 

r I z 



(23) 



(24) 



is an eigenfunction of A s with eigenvalue —1(1+1) = —2. 
In ([TBI we will choose so that 



2 *(X) - Ba{xd y - j/9 x )*(X) 



(25) 



is a multiple of *F Let 



n 



( 


-CJl 




























-U> 2 


I 











Substituting into 

f! 2 *(X) - 5ct(.tg\, - yd x )*(X) 



(26) 



(27) 



produces 



B)x 



(28) 



V 



o 



This will be a multiple of ^ whenever u>2 
So 



f -co 

w 

_0_ 

V 



— VcJ 



I? 



(29) 



/ 



We put this form of Q into the wave equation to get 
an ordinary differential equation for w: 



—w' j w + d{ w ) + i^ 2 ~ B)w 



(30) 



This is known to have exponentially localized solutions 
provided lim s ^ g(s)+( f- B)s < (i). So we require 

— p 2 + uj 2 — B < (where lim s _>o = wmcn i s tnc 
case for all lo with < to 1 < p 2 + B. Then, according 
to the results of 0], there are C 2 exponentially localized 
solutions to (|30p . The corresponding functions u(X,t) = 
e ta, &(X)w(r) are solutions to Q. 

Using the form of Q we now show that a solution of the 
form u(X , t) = e ^f(X)w(r) has nonzero spin parallel to 
the direction of the magnetic field B. The component of 
the spin parallel to B is 



S z [u] 



w 2 nii(X) ■ (xd y - yd x )^(X)d 3 X{31) 



> 



(32) 
(33) 



if lo > (and both x and y are nonzero). So, 



SJu] 



^(x 2 + y 2 )d 3 X 



(34) 



is nonzero. 
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